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ABSTRACT 

In a recent paper* by T. J. Harvey, there is presented the idea of 
employing known solutions of free vibrations of nonlinear single-degree 
conservative systems to obtain solutions of forced vibrations of these 
Systems. The forcing function is taken to be of the same form as the 
free response and the solutions are obtained quite easily. Harvey makes 
the point that while this excitation is not the sinusoidal excitation 
that many analysts have attempted to study, for many purposes there is 
reason to prefer some alternate form of excitation if by doing so useful 
solutions are readily obtained. Ina paper of several years back, M. A. 
Sadowsky* derived the characteristics of several nonlinear spring systems 
for which the free vibrations could be obtained in rather simple closed 
forms. In particular, Sadowsky introduced two new springs which give solu- 
tions in terms of the elementary functions and several others which give 
solutions in terms of the Jacobi elliptic functions. 

Since there are actually only a few fully developed cases of forced 
nonlinear vibrations available to analysts for testing the power and ac- 
curacy of analytical methods and for other purposes, it appeared to be of 
interest to apply the viewpoint of Harvey to several of the particular 
cases given by Sadowsky so as to work out all the details of the corres- 
ponding cases of forced vibration and thus afford several new examples. 
The major content of this thesis is the full development of these cases. 

In addition, however, the viewpoint of Harvey is applied to two cases 


of piecewise linear springs simply for the purpose of illustrating such 


*See references on page 4l., 
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applications and for the purpose, as well, of obtaining additional cases 
of nonlinear vibrations. 

Also, since the forcing functions differ from the sinusoidal func- 
tions that many analysts have used, it was considered desirable to obtain 
and exhibit a simple measure of the degree of "distortion'’ of the actual 
forcing function as compared to the sinusoidal shape. Likewise, the dis- 
placements are not sinusoidal functions of time, and the "distortion" of 
the response (as compared to a sine curve) is also examined. 

This work was done largely during the spring of 1961 at the U. S. 
Naval Postgraduate School, Monterey, California. The writer wishes to 
express his appreciation to Professor J. E. Brock for having suggested 


the topic and for guidance during its development. 
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SYMBOLS 
time = independent variable, chosen so that (f/0) = 0 
dimensionless independent variables 
deflection, normalized 
external force applied to the system 


spring function such thatS/o)= 0 and $0) = 1.0 
S(#) has the same dimension as does 4 


assignable dimensionless constant as defined in the relation (5) 
mass or moment of inertia of the system 
spring constant for a small deflection 
Ky, = natural frequency for a small free oscillation 
—— = period of a small free oscillation 
period of a free vibration Ge = 0) 
frequency of vibrations 
ak = period of forced vibrations 
4(4) = amplitude of deflection 
denotes the evaluation at f= A 


modulus of Jacobi elliptic function when a solution contains 
such functions (refer to Table 1) 


modulus of Jacobi elliptic functions wher a spring function 
contains such functions (refer to Table 1) 


the complete elliptic integral of the first kind with modulus % 
the complete elliptic integral of the first kind with modulus ~y, 
the distortion of #) wave form as defined in Eq. (17) 
the distortion of Fit) wave form as defined in Eq. (18) 


a constant for the springs of cubic parabola types such as 


RY) =Y+p VY? 


deflection which is not normalized 
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Ll. Introduction. 

The equation We + ht = Eevee (1) 
describes the behavior of a linear single-degree-of-freedom conservative 
system subjected to a sinusoidal excitation. 

A natural extension for the study of a nonlinear conservative system 
is 

m4 + kS(t) =F an wt (2) 
where S|¥) is a nonlinear spring function, such as Six) = 444° . Many 
investigators have studied Eq. (2) in detail and obtained solutions which 
indicate the effects of parameter {%) upon “A . But Eq. (2) is complicated 
to solve analytically, and in most cases the solutions are obtained by 
approximation methods. 

T. J. Harvey (ref. 1) suggested a simple analytical method of study- 
ing the forced vibrations of nonlinear systems by making use of the solu- 
tions for free vibrations, i.e., where the right hand side of Eq. (2) is 
zero. Some years ago M. A. Sadowsky (ref. 2) presented simple solutions 
of free vibrations for various nonlinear springs. 

This thesis applies Harvey's method to several of the cases given by 
Sadowsky, thus obtaining several interesting solutions for forced vibra- 
tion of nonlinear systems, 

In addition, a study is done on the following: 

(a) The wave form of (f(t) is compared with sinusoidal one which 
is that of a linear system. 

(b) Traditionally the right hand side of Eq. (2) is a sinusoidal 
function, as shown, By Harvey's method it differs from a sinusoidal func- 
tion, but is still a periodic function. These differences are compared. 
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(c) Harvey's method can be applied to other than Sadowsky's 
nonlinear springs. It is applied to the spring characteristics which 


consist of linear segments, as examples. 





hesic ecuaétions, 
The forced vibration of an undamped nonlinear single-degree-of- free- 
dom system is described by the equation 


wig + f Sf) = FM) (3) 


4 


where S/») is a nonlinear spring function normalized so that Sicj/=o 


1¢ 


and Sjoj= 3.0 , and Fr) is the external force applied to the system. 


4 


Eq. (3) is simplified by substituting: 


oo ky, and L= &ft (4) 
en a Se mn F () a) 


Eq. (3a) as it stands, is of great difficulty. Even if F(t) is 
periodic and if we presume (by way of ana logy with linear systems) that 
slight damping present in the system but not accounted for in the mathe- 
matical description causes the "transient" to die out, the "steady state" 
Solution is generally not obtainable except by difficult analyses leading 
to approximations. However, the analysis is greatly simplified by henceforth 


confining attention to only such f(t) as satisfy the relation 
F(2) 
~ & SL HO] 


where P is a dimensionless constant, either positive or negative. Such 


? (5) 
a forcing function is called a "natural forcing function" (ref. Ll). A 
positive P means the external force is in-phase with restoring force, a 
negative ‘e means out-of-phase (180), and P=C means free vibration. 
Substituting the relation (5) into Eq. (3a), we obtain 
dy 
7S + (im p) Slt) = 0 (6) 
This is further simplified by putting 


g == fi-/ Tc (7) 








» —_— - 
en nee pa Lt (8 


which is described as the case of the pseudo-free oscillation. To solve 
Eq. (8), let the initial conditions be 
fic) =o, =o = Ve 
The potential energy stored in the spring is 
Ulf) = J Stdp ue) =e 
Furthermore, by conservation of energy, Vo is related as 


Ve = 2 UF) 


where 9 is the amplitude of deflection. If T/4 is defined as the quarter- 


o£, 


period of oscillation, 4 can be writtenas ¢= 7( Vj.) . 


Eq. (8) is integrated and simplified as 


- 


(ZY = 2u 9-24) ) 


Integration of Eq. (9) leads to 
¢ dy ‘ 
C= } Pup 2ugpye sis 
Again Eq. (10) can be written in the following general form with 4 
as a parameter, 


6 = O(¢. #) (11) 
The inverse form of Eq. (11) can be expressed as 
¢ = ¢(6,#) 
Since the dimensionless variable © is closely related to time t, 
using Eqs. (4) and (7), 
@ =4,t[I- Pp (12) 
Then Yom f(a t lip i ) (13) 
This is the solution of the forced vibration of Eq. (3) in which the 
external forces F(t) satisfy the relation (5). For free vibrations /(t)=< 


and Ome» Then the solutions of free vibrations are from Eq. (13), 





Comparing Eqs. (13) and (14), the solutions of forced vibrations are 
obtained from those of free vibrations by multiplying the independent vari- 
able t(= time) by the factor Vi-=p - The assignable dimensionless con- 


stant . is the ratio of the external force to the restoring force, namely 


isis 
R S(#) 

where the relation (5) is rewritten. The amplitude of external force F 
is given as 2.0 k or 3.5 , for example, which means the force unit corres- 
ponding to the value of the spring constant (at small deflection) multiplied 
by 2.0 or 3.5 units of deflection respectively. Or in other terms, Eq. (3a) 
represents a forced vibration wherein S(t) plays the role of a restoring 
"force" and = F(c) plays the role of an exciting force; the time scale 
has been chosen so that both ''forces'" have the same dimension as the "de- 
flection response", 7(t). Thus we can say that the amplitude of exciting 
force ( ER ) has magnitude 2.0 or 3.5, say. For every amplitude of de- 
flection Ci , the values of 5S = S{¥) are found, which, in turn, deter- 
mines the assignable constant p ; 

Sadowsky's paper gives the essential details for ten cases of nonlinear 
springs for which the free oscillations are relatively simple. Two of the 
solutions are in terms of elementary functions, and the others involve the 
Jacobi elliptic functions. Corresponding to each of these cases, by using 
Harvey's method, it is possible to construct a complete solution for the 
forced oscillations resulting from "natural" excitation. In this thesis 
four of the cases are developed in detail--the two cases involving only 
elementary functions and two of the cases (cubically hardening and cubi- 
cally softening -- i.e., Duffing's spring) involving elliptic functions. 
Also, two cases of piecewise-linear springs (cases not included among those 


given by Sadowsky) are also developed. 
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Sadowsky's solutions are expressed in terms of the initial velocity 
Ve , which can be replaced by 4 ,» for convenience, with the energy rel: - 
tion \.= 2U(#) in the conservative systems. The solutions are tabu- 
lated in Table 1 and their approximate spring characteristics are shown 
in Fig. 1. The symbol O& represents the modulus of Jacobi functions, K 
means the complete elliptic integral of the first kind with modulus ~™ 
In Fig. 1, the numbers in circles represent the spring numbers as in- 
dicated in Table, Springs marked as 1, 2, 3, and 4 are drawn heavily. 
The family springs 7 is in between springs 1 and 4 depending upon the modu- 
lus A, , where for %=0O it coincides with spring 1 and for %,=|-0 
it coincides with spring 4. The space between springs 4 and 3 is filled 
with the family of springs 6 depending upon the modulus ™% . For o%=0 
it coincides with spring 3, and for %i=/.0 it coincides with spring 4. 
Each one reaches a maximum and drops down to zero. The family of springs 
5 is of softening type, and so is the elementary spring 2 obtainable with 
A =1.0 . Springs 8, 9, and 10 are cubic parabolas, they are not shown 


in Fig. 1 since they can readily be visualized. 





Table 1 - The solutions of forced vibrations for Sadowsky's springs 
(adapted from Sadowsky; ref. 2) 


The solutions of forced vibrations such as 
4 : : , SN ate —_—— ic 
Ji{ at + Fo { a aa F(c) 


are tabulated below, in which Etc) satisfies the relation (5). 
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Fig. 1 Characteristics of Sadowsky's springs, 
(from Sadowsky, ref. 2). 
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Fig. 2 Distortion of the Fig. 3 Distortion of the 
wave form. A is wave form. B is 
defined in Eq. (17). defined in Eq. (18). 





A. Frequency response, 

When the exciting force is such that the relation (5) is satisfied, 
the period of forced vibration is readily obtained from the solution of 
free vibration. If the period of free vibration is obtained as: 

T. = TOY) 


then by comparing Eqs. (13) and (14), we see that 


2 = ae 
Sea (15) 


ni = mal 
Since /o is varying with 4 , the ratio baal may be related to the 


period of a small free oscillation _— 





T -(2) —L 
TH Ce Vt-p (16) 


This means that corresponding to a single frequency response curve 
on Ff versus (4) , obtained for free vibrations, we obtain a full 
n 


family of response curves for forced vibrations, i.e., curves of f versus 


(=) for various values of fr : 





B. Wave form of /7/(t) . 

When the exciting force is such that the relation (5) is satisfied, 
the solution is exactly the same as that for free vibration except for 
the change in frequency, which is JI-/ times that of the free vibra- 
tion. Therefore the hae form of ¥(z) depends upon the amplitude of 
deflection ¥ only, and not upon the amplitude of exciting force. When 
4 is small, the springs are very nearly linear, and the wave forms 
are almost purely sinusoidal. As ¥ increases, the wave forms of th) 
are distorted from the sinusoidal forms. 

Furthermore when the spring characteristics are symmetrical, the 
wave formS are quarter-wave symmetrical. In addition, they tend to lie be- 


tween a rectangular and a triangular wave form, 


To show the distortion of Ga) wave forms, a measure A is defined as 


ee a 
A = 2(2)/ & dt (17) 


If we let the quarter period T/4 correspond to % , the measure A 


represents an area under the curve of & as shown in Fig. 2. Note 


that A=], for a linear spring or a small oscillation. 
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C. Wave forms of F a) ; 

When the external force F(t) satisfied the relation (5), the wave 
form of F(t) is the same as that of the spring force S| KM] . Therefore 
it depends on 4 only, not the amplitude of the external force F ° 
For a symmetrical spring, the wave form of Fit) ts quarter-wave symmetri- 
cal in the same manner as that of f(t). For a small oscillation [ (7) 


is almost a pure’ sinusoidal form, from which the wave form is distorted 


as ¢ increases. 


To compare this wave form with the traditional forcing function of 


sinusoidal form, a measure B is defined as 


B=24@) | #4 


(18) 


Note that: Bae5 for a linear spring or a small oscillation. 


RE 





3. Examples. 
The results of six examples are presented in Fig. (4) through Fig. (15). 
Figures 4, 6, 8, 10, 12, and 14 are all similar representations as follows: 

(a) Ordinates are jf (or 4h where € is a constant). 

(b) Solid curves represent frequency reSponse curves for 
labled values of [ = F(%) » including F=0 , the case 
of free vibration. 

In-phase forced vibrations are indicated by positive values of F , 
out-of-phase by negative values. These curves relate 4 (ordinate) to 
4, (bottom scale of abscissa). 

(c) Distortion of ¢(t) is indicated by a measure A, distortion 
of F(t) by a measure B, both functions of 4 alone, and 
plotted on vertical scales. 

(d) The dotted curves show spring characteristics SF) (top 
scale of abscissas) as functions of 4 eine bags eeiee 
and 14, this relation is shown to smaller scale in the 
inset and not in full scale as in the preceding figures. 

Figures 5, 7, 9, 11, 13, and 15 actually display the distorted wave 


forms for ¥t) and F(t) . 
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A. Example 1 j= =~ 
The characteristic of the spring is extremely hardening as shown in 
Fig. 4 (dotted curve). The solution of the forced vibration is already 


tabulated in Table l, namely 


_ 


MH Y oe att $ dan Wl (19) 


where 


WY =(Wy sec E) TP, p = JE (20) 


(a) Frequency response, 


From Eq. (20) we can get the solution which indicates the 


7 27 : 
effects of frequency upon 4 , making use of |= ay and h= a ao elit 


the dimensionless form, we get 


Sere. Ce a 
ae Fra (21) 


For the given value of F , / is obtained by dividing 
F by S(#) , then the ratio 14. is obtained from Eq. (21). The fre- 
quency response curves are plotted in Fig. 4. 
(b) Wave form of f(t) . 


Eq. (19) is rewritten in a convenient forn, 


: ain’ (tind pin wt) — 
s = ene ane) (22) 
For a small oscillation, Lem, — = An Wt 


With Eq. (22) the wave forms of #{t) are plotted for 4 = 0, 
0.5, 1.0, 1.30, and % in Fig. 5(a). Since Eq. (22) shows that 7(t) is 
quarter-wave symmetrical, plotting is done for a quarter period. For a 
small oscillation, it is clearly a pure sinusoidal form, and as a in- 
creases it is distorted. To show the distortion, the measure, A{as defined 
in Eq. 17) is found. The integral is not easily performed analytically 


so that a numerical integration method was used. The measure A is a 


I3 





—se 


function of 4% , and is plotted along ¥ in Fig. 4. 
(c) Wave form of //t) . 
Using Eq. (19) and the relation (5), 
=F lee ee tc = Bw 
Fae oR S[ 4) and F = PkSCF 


rearranging these into a dimensionless form, we get 








Fit) aenwt 
ee eg (23) 
j Flt ? : 
For a small oscillation, o = dint 


which is purely sinusoidal. As 4 increases, it is distorted. With Eq. 
(23) the wave form of F(z) is plotted in Fig. 5(b). 

To compare F(t} wave form with the traditional sinusoidal function, 
the measure B(as defined in Eq. 18) is obtained: 

ay, 
B= E()[" 20 tt = cea7 
which is a function of # only, and is plotted in Fig. 4. 
(d) Discussion. 

This spring has an extremely hardening characteristic, At 
¥o_ , the spring force becomes infinity, its slope is infinity, and 
the measure B is zero. For slightly smaller values of 4 , the messure 
B is far away from unity. Fig. 5(b) shows that the wave form of F(t) may 
become quite distorted. 

When A is far away from the value of 7 , the measure B is not far 
away from unity, and Fig. 5(b) shows the wave form of /(t) is not so much 


distorted, 
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B, Example 2. oa ~ — 
This spring has a softening characteristic as shown in Fig. 6, 
The solution of the forced vibration is already tabulated in Table 1. 
dah $ = denh % ain uit (24) 


where 


w= (a ack F fio ee se (25) 





RS 
(a) Frequency response. 
ida eae 
From Eq. (25) ee = coh ja (26) 


For a given value of 5 ; e is obtained as a function of 
ve , and the dimensionless values of are found. The frequency re- 
sponse curves are plotted in Fig. 6 forfF=e and tof. 

(b) Wave for of f(t). 


Eq. (24) is rewritten in a convenient forn, 
git) _ _denk\( aich F ain wt) 
$ g 


? . ; - 
For small oscillation, hn ae) = Un WIG 


(27) 


which is purely sinusoidal. With Eq. (27) the wave form of tt) is 
plotted in Fig. 7(a) for ta 0 7 1.0, e2n0, 5.0, and OO, As 4 increases 
it becomes distorted. And as i approaches infinity it goes to a rec- 
tangular form. 

To show the distortion, the measure A is found and plotted along # 
in Fig. 6. The integration to find the measure A is performed by a 
numerical method. 

(c) Wave form of F(?#). 
Using Eq. (24) and the relation (5), we get 
F(t) = pS,  F= pk S08) 


rearranging these into a dimensionless form, 


A) oe 
FT (1 ~ Leh?eoctat Ye 
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For a small oscillation, St me dn ol 
which is purely sinusoidal. As if increases it is distorted, With Eq. 
(28) the wave form is plotted for ¥# = 0, 0.5, 0.66, and 1.0 in Fig. 7(b). 

To compare F(t) wave form with a sinusoidal one, the measure B 
is found, 

B= cwsh F (29) 

which is plotted along 4 in Fig. 6. 

(d) Discussion. 

The spring is of softening type: at # = 0.66 S reaches 
its maximum value (which is equal to 0.384). Beyond this point it decreases 
and the spring force approaches zero as if approaches infinity as shown in 
Fig, 6. Up to + = 0.66, the measure B is not far away from unity, and 
Fig. 7(b) shows that the wave form of F(t) is not distorted appreciably. 

Far beyond the point 4% = 0.66, the value of B is much greater than 


unity and the wave form is greatly distorted. 
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C. Example 3. Duffing's equation. S(#) — g + =f 

When the spring characteristic is a cubic parabola, the equation of 
forced vibration is 

mY + ke Y+sY") = F(t (30) 

in which Eq. (3) is rewritten here. Y is the deflection. If it is ex- 
pressed in inches, p is in incre and k in lbs/in. 

To solve Eq. (30) with the right-hand side set equal to zero, the 
general cubic parabola Oey sy can be expressed in the form of Table l, 


ie, y (f+27) by proper substitution regardless of ff . 





ee. 
Let Caer | (31) 
Then the equation of free vibration becomes 

wa ll — 

4 +Un(G+2P¢ } 


in which ¢ is a dimensionless quantity. 


Let S(¢) = rs at (22) 
Then we obtain the solution of forced vibration immediately from 
Table 1. 
Wt) = % On(ad-K , & ) (33) 
where W = w,(Its DIP 
x = Cre as the modulus of Jacobi functions (34) 
P= Fs ie in a dimensionless form 


And K is the complete elliptic integral of the first kind with modulus 


% . 
OX » Le., K ae { (1- ang)? de 


(a) Frequency response 


From ~~ (33) and (34) with T= » we get 
ae ee (Ze 
ER = 


In Fig. 8 the frequency response curves are plotted for F=0 ok Be : 


and +2 $f 
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(b) Wave form of Hae 


Eq. (33) is rewritten: 


Att} —s | 
So = Cy (4tT-K,X) 
f . : era 
For a small oscillation, a = dnl 


which is purely sinusoidal. As a increases it becomes distorted. 
However there is a limit in the distortion. As ¥ approaches infinity the 
value of & does not exceed a - In Fig. 9(a) the wave form of (f(t) 
is plotted for PI = Q and # = 03, 
To show the distortion the measure A is found, 4 =(.) ew 

which is plotted in Fig. 8 along 4 » It does not change appreciably 
from unity. Even if 4 approaches infinity, A approaches 0.942. In other 
words, the wave form of 4{t) is not distorted much from the pure sinusoidal 
form, 

(c) Wave form of F(t) « 


From Eq. (33) and the relation (5), 


Fin =p S[4] FF = prs SCH) 


rearranging these into a dimensionless form, we get 


Ee. | n(wt-K) + tf (ut-K) | (36) 


F [+i 7 ; 


with the medulus X 


ee CS, 


— 


For a small oscillation, goo fF = AMwt which is purely sinu- 
soidal, With Eq. (36) the wave form of F(t) is plotted for 4 = 0, “tec. 
2.0, and infinity in Fig. 9(b). 


To compare [-(t) with a sinusoidal force, the measure B is found, 
= t 
7 x (i 2 
B — 2K 29°44 
which is plotted along 4 in Fig. 8. 
(d) Discussion. 


This spring is a hardening type as shown in Fig. 8 (dotted 


curve). Every type of cubic parabola, such as Y+sy’, can be reduced to 


22 





the form S(t)= f+ 3 a by letting ferme: without losing the 
generality of the governing equation (Eq. 3). 

Fig. 9(b) shows that the wave form of F(t) is not distorted 
to an appreciable amount, and the measure B is not far away from unity. 

The application of Harvey's method to Duffing's equation 
was considered in a recent note (ref. 3) which appeared after this thesis 
was nearly completed. For the hardening spring, considered in this section, 
the results given in the note agree with those given here. However, in the 
case of the softening spring, treated in the next section following, the re- 
sults do not agree since the note contains an error. In place of the cn 
function in the note, the cd function should appear for the softening case, 


and the argument of the function should also be slightly modified. 


Ja 


diel 


> <4 





= See - 











D. Example 4, Duffing‘s equation. S (4) = 4-34 
When the spring characteristic is a softening type of cubic parabola, 
the general characteristic equation is expressed as 
Ri) = Y- pY 
which can be reduced the form of Table 1 by letting Y==s? , where 
¢ is a dimensionless quantity. Then without losing the generality of 
the governing equation (Eq. 3), we get 
S(#) = 4-3 
The solution of forced vibration is already tabulated in Table l, 
it) = & An(wt, x) (37) 
where 
=e F) HP 
of 
A= 94-7 (38) 
p = 45 ap 
(a) Frequency response. 
From Eq. (37) and (38) with T= 4K we get 
Te OE) (ES) on 
where K is the complete elliptic integral of the first kind. The fre- 
quency response curves are shown in Fig. 10 for fF = 0 and +02 Mes ° 
(b) Wave formof ¢/t). 
Eq. (37) is rewritten in a convenient form, 
ot A (uit, x ) where & =%(4-7)% 
For a small oscillation, en, at = men wt 
which is purely sinusoidal. As 4 increases it is distorted as 
shown in Fig. 1ll(a). 


To show the distortion the measure A is found, 
a” 
IL 2X a ig 
A as eza, by ‘cr 1— (40) 
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and ig plotted in Fig. 10 along 4 . 
(c} Wave form of F(t) 
From Eq. (37) and the relation (5), 


= An. = 9 eee 
Fit) = pre S[H?] ana F =P = SF) 


ry, 





rearranging these into a dimensionless form, we get 


F(t) { na ae 4 
Te pono | a(wie «) - ze An (Wt, | 
z 2 ¢ (41) 
11 oscillati hin KO. ainutt 
For a small oscillation, Goo | 


aed 


which is sinusoidal. As g increases it is distorted as shown in Fig. 


110d). To compare this wave form with a sinusoidal one, the mea- 


Sure B is found, 
= ne eo ae : 
B =k) Ait Ry" + (4) (A%) ['- aol a | (42) 
which is plotted along 4 in Fig. 10. 
(d) Discussion. 

This spring is of softening type, and at $= G.82 it reaches 
its maximum. After this point it drops down, and at Yalz the spring 
force becomes zero. 

The measure B does not change appreciably up to #= 0.82, 
and Fig. 1li(b) shows that the wave form is not distorted much. 

For any value of f# which indicates the nonlinearity of the 


spring, the above solution can be applied. 
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E. Examplq@ 5. 
Suppose 4 symmetrical spring 
(as shown) is described by the equations 


lal<c 
eas 


Sip) = hig 
St) = Rep +( ker ka) 


and it is considered the case of k)Oonly. 


The equation of free vibration is 


mie + S(t) = © 


Piecewise-linear softening 


3 


Spring®, 





(43) 


When the amplitude of deflection # is less thanc , it is linear, 


and solution will not be included here. 


ed, 
Let us denote #(t) for 
4 (D) for 
where Zc = the time when Y(t.) = Cc 
and let (= by and ur = ken 
Then Eq. (43) becomes 
$+ Wy 9 
+ ih t(Wnasje = 
The general solutions of Eqs. (43a) 


¥(t)= D An(ytt p) 


f(t) = Dz wen (wat +b) + (I- Fe) 


The case of §>C will be consider- 


oO (43a) 


are: 


LETSYG (44) 


* In section E and F (which follows), for convenience in notation, 


only the first quarter cycle of motion is considered, 
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where D, DEE P » and @ , are arbitrary constants, which together with 
Te may be determined by the following boundary conditions; 
4(c) = O as an initial condition. 


Wo) = ez u( #) by the conservation of energy. 
f(t) = f(t) =¢ 
¥(t) =£(t) 


Solving for the constants in Eqs. (44), we get 
a i ee ae 2 ‘2 
rp) = [ c +2¢e(¢-¢c) + 5.(F-c) | 
~ See ee - 
D, = -€C|- [ee (45) 
te = ~ ten = 
| Ue) _ ht. $ =0 
d, = Hin” Gs) - 7 ee 
To obtain the frequency response curves, the period of free oscil- 
lation ( p = 0) J, can be found: 
t=, £X)q 
because fall) reaches its maximum value of 4 at az. o solving for J, 
7 : ee ey) 
we get: La ZB (= - g, ) 
_s ) (1 sz 4) 
=G; 1, 
where aM the period of a small free oscillation. 
By making use of Eq. (16) for forced vibrations, the period of forced 


vibration is obtained, 





EZ = (%)('- be) = (46) 
As a specific example let us take b=4h, in the above symmetrical 
spring. Then we get (% =2wW, , and the solutions of Eqs. (44) become 
Gt) = D, amust alates 
42(t) =(#+3c) dn (azt+,) -3¢ fo eA 


Dearie + 2am fc) te(F-c)] ts 


(44a) 


where 


| 





Ye ee 
Co = yy AM ( a 


i ea 
Fz) mee Ken (Fz) 


(a) Frequency response. 





A, 
re - = = putt 


From Eq. (46), 
ee 
Txt aS 


where & is an assignable constant as defined in the relation (5). [In 





(468) 


Fig. 12 the frequency response curves are shown for fF=0 , +t 0.5 8,0 


and + ke. 


) 


(b) Wave form of A(t) : 


Eq. (44a) is rewritten. 


ith «2 man wT oft ST. 
j Bc ~~ s 
site = (144 35) ain (WT +4.) — FZ L675 


For ¥<C the Sam is purely ei da, and as ¥ in- 
creases it is distorted, At around ¥#=2 5c , distortion reaches a maxi- 
mum, and for further increase in 4% the wave form approaches a sinu- 
soidal form again. In Fig. 13(a) these two wave forms are shown: one is 
Sinusoidal and the other shows the maximum distortion at 4 =2-5¢ : 


To show the distortion the measure A is found, 





A= B (1 coouste ) gg Zt (1+ Apeeleth rd.) 
- EFI (I- sh 


which is plotted in Fig. 12 along tfc 


(47) 


- Actually the distortion is so 
small that the measure A is not much greater than unity, where the maxi- 
mum measure A = 1.023 at 4=2Se : 
(c) Wave form of F(t) . 
From Eq. (45) and the relation (5), and with the spring 


force -S(¢) = (3c +4) b, l#|2c , for this spring, we get 


Be 





F(t) 4D, 





a roe din We f CL TST. 
ee ; Di a a 
Eo eee StS % (48) 


| 
In Fig. 13(b) the wave form is plotted for ye ; 4% =2C 


and $c . Since the spring force has a slope discontinuity at eae, 
the wave form of F(t) is slope discontinuous at t=f.. 

For fac it is a pure sinusoidal form. As ¥ increases 
it becomes distorted; maximum distortion occurs at p=? fe approximately. 
Further increases in if bring out a sinusoidal form again, similar to 
f(t) wave forms. Every wave form of f(t) does not go beyond a limiting 
curve which is marked as "Envelope" in Fig. 13(b). 

To compare the wave form of F(t) with a sinusoidal forn, 


the measure B is found, 





4D >, eee VI icone 
6 ~ F430 (! coat. ) % -@, t Y, - Cod ( Ul. re.) (49) 


which is plotted along UVA in Fig ..(l2- sac be =2-5¢ , approximately, § 
reaches its maximum value of 1.22, At this point the distortion is a 
maximum. 

(d) Discussion. 

The solution of free vibration is found by a classical 
method, and for the forced vibration Harvey's method is applied to that 
of free vibration. 

Fig. 13(b) shows that F(t) is not far away from a pure 


Sinusoidal form, and the measure B is not much greater than unity. 
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F, Example 6. Piecewise linear hardening spring. 
Let us consider a symmetrical 
spring as shown, Let Ak, ‘a k, then 24<-4, 
For free vibrations, 
m 4 + S(f#) =O 


Let the initial conditions be 





d (c) =O ; ¥ (0) — 2 ulF) 


Then the solution is obtained from Eqs. (44) and (45), 


40 =D dn vat 94TSk 
= a | a. (50) 
4(t) =(F- 48) Un aT h,) + £6 Late 
where . , ae 
D =[c +2e(G-<) Fal F<) 
! —- Cc 
Fe = | Ain"(S) 
= Ken Coe ee Sm ESS 
ae Mlk Ge) 2AM Gs 
(a) Frequency response. 
From Eq. (46) and with 2”# =e. , we get 
Ty... oe 
a > ) a (51) 


In Fig. 14 the frequency response curves are piotted for 
F #0, ¢ hc , and +5hC. For gee , 1t is linear and the ratio 
( ivr, ) is unity. And as ib becomes much greater than C , the spring 
characteristic again approaches linearity and the ratio ( ty, ) approaches 
0.5. 
(b) Wave form of F(t) . 


Eqs. (50) are rewritten: 


dD 
a = ee Ain iE Es eas (52) 
Sol ; 
Ss = (is <3) un ( wt + @, ) ee Lats Y, 


in Fig. 15(a) the wave form of g@) is shown for ¥Soc 


and esiSl. For £<C , the spring is linear and the wave form is purely 
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Sinusoidal. When 4 becomes greater than C , the wave form ig distortei, 
at approximately & = he the distortion reaches a maximum, For further 
increase in ¥ , the spring characteristic approaches linearity again, 

the wave form approaches a sinusoidal form. 


To show the distortion the measure A is found, 


T ay 1, Fe = ; ) # 
A= rE GY (Ie ewt) + ay rds) esate) 
op 50 \ ol ae pers 
+ 2(25)(i- Hea) oT 


which is plotted along th in Fig. 14, At =i sc the measure A reaches 
a minimum (A = 0.965). 
(c) Wave form of F(t) . 
The spring characteristic is 
S{F) = hg - Phe for |f|2- 


Substituting Eqs. (50) into the above equation, we get 








fo = - den ait e<etet. 

aa “+ # -3C 

Fa). gin(at+ %,) tests % (54) 
E 


In Fig. 15(b) the wave form of F(t) is shown for ¥<cC_ , 
Pi = HH , and 3.0C ., There is a limit, beyond which the wave ferm 
does rot go, as indicated as "Envelope" in the figure. 


To compare the f(t) wave form with a sinusoidal form the 


measure B is found, 





D, : h } i V/2 : - | 
f3 a qaoe (| hee) we, UE caer Cod (4a ke r@, ) (55) 


which is plotted in Fig. 14 along #4. For #<c the spring is linear and 
F(t) is a pure sinusoidal form, so the measure B is unity. As ¥ ZOES 
beyond c , the measure B decreases rapidly, and at g=lSc, approximately, 
B = 0.735 as a minimum. For further increase in 4 , the measure B 


slowly approaches unity where the spring characteristic becomes linear again, 
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(d) Discussion. 


This example is similar to the previous one, from which 


the general solution is obtained. The measure B and the wave form of 
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